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Abstract 

We investigate the critical parameters of an order-disorder quantum phase 
transitions in the spin-^ J —J' Heisenberg and XY antiferromagnets on square 
lattice. Basing on the excitation gaps calculated by exact diagonalization tech- 
nique for systems up to 32 spins and finite-size scaling analysis we estimate 
the critical couplings and exponents of the correlation length for both mod- 
els. Our analysis confirms the universal critical behavior of these quantum 
phase transitions: They belong to 3D 0(3) and 3D 0(2) universality classes, 
respectively. 
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The equivalence of the critical behavior of D dimensional quantum spin systems and 
D + 1 dimensional classical spin systems is well recognized. This idea, combined with finite- 
size scaling was used previously many times to discuss the critical properties of infinite 
spin systems, see e.g., and for a review. However, these investigations were strongly 
limited with respect to the system size. Due to the recent advances in computer technology it 
is possible to treat bigger systems, e.g., up to 36 sites for spin |, and consequently to extract 
their critical properties using finite-size scaling method Our aim is to present in this 

paper the results of such an investigation of critical parameters (coupling and exponents 
of correlation length) for the J — J' Heisenberg and XY spin-^ antiferromagnets on square 
lattice. 

The Hamiltonian of the model whose critical behavior is examined is given by 

H = J s,-s, + j' Yl Sk-Si. (1) 

<i,j> <k,l> 

The first sum, denoted by < i, j >, runs over pairs of nearest-neighbors on the square lattice 
connected by thin bonds (see Fig.l), whereas the second one, denoted by < fc, Z > — over 
nearest-neighbors connected by thick bonds. In the case of Heisenberg spin system three 
Pauli matrices are included into scalar product in Eq. (|I]), in the case of XY system — 
only two. The model represents an antiferromagnet, i.e., both couplings are positive and 
additionally, J' > J. Clearly, what one can see here is the competition between long-range 
Neel order and the tendency to the formation of local singlets of two neighboring spins, 
coupled via J'. In limiting cases this model reduces to the long-range ordered Heisenberg 
antiferromagnet on square lattice for J = J' on one hand and on the other — to disjoint 
singlets (no staggered magnetization) for J'/ J ^ oo. At some (J'/ J) there exists a quantum 
phase transition between the gapless Neel phase and a gapped 'singlet' phase (quantum 
paramagnet). The properties of J — J' Heisenberg model on the square lattice were first 
examined by series expansion (SE) p| and more recently by renormalized spin wave (RSW) 
approach exact diagonalization (ED) and coupled cluster method (CCM) ^ in order 
to observe the interplay between the local singlet formation tendency and the long-range 
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Neel order. Although, in general, all those methods predict the existence of quantum phase 
transition, they differ drastically in the estimation of the critical coupling. Furthermore, 
there was only one attempt |^ to find critical exponents, however the error of this estimation 
was rather large. 

Let us now rewrite the Hamiltonian in such a form that the two above tendencies will 
be seen more clearly: 

H = g{ ^.-^. + 4 E 'Sk-^i), (2) 

<i,j> ^ <k,l> 

< g < 1. Note that in Eq. (0) one has the same summation specification as in Eq. (||). The 
coupling constant g determines the relevant energy scale in the model under consideration 
and the term l/g"^ = X is analogous to an inverse temperature. 

In what follows the estimation of the critical value of gc and critical exponent l/u for 
the Heisenberg and XY Hamiltonians will be described. At the beginning one finds by the 
exact diagonalization (Lanczos algorithm) the spin-gap, defined as A = — i^o (-£'0(i) is the 
lowest energy in the 5*^ = 0(1) sector) dependence vs. X = l/g"^ for the Hamiltonian (0) on 
a square lattice for a sequence of finite systems. Note that all the systems, being elements 
of the sequence should be invariant under the same symmetry operations; in the opposite 
case it is not possible to find a proper scaling. In the system under examination one has 
only three such systems: with = 8, 18, 32 spins. Two of them are shown in Fig. 1, the 
third one has the same shape. 

The spin-gap, multiplied by the linear dimension of the system, A^/N, allows one to find 
the pseudo-critical points 0]: for the Heisenberg and XY systems one has three such points, 
collected in Table 1 (Heisenberg) and 2 (XY). Next one should extrapolate the sequence of 
pseudo-critical points to infinity. However, it is not possible to employ here any algorithm 
improving the convergence of the finite- lattice data sequence the series is extremely short. 
Therefore we find it more accurate to use graphical methods to find the critical value of 
g. The sequences of pseudo-critical points for Heisenberg and XY system plotted vs. 
are shown in Fig. 2. The estimate of l/g"^ = 2.46(2) for the infinite Heisenberg system 
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should be compared to the value 2.56 obtained by SE p, to the value of 3.16 obtained by 
CCM method, to the value of 2.45 from ED [|] and finally to the value of 5.0 from RSW 
approach. The extrapolated value oil/ g1 = 4.56(2) for anisotropic XY system is higher, as 
one should expect, than that for Heisenberg system since anisotropy act aginst the singlet 
formation |Q. 

The critical exponent for the correlation length, z/, may be estimated from the behavior 
of the Callan-Symanzik /3-f unction [ll0|,|3|,[l[| 



/^(A)A=^lnbA(A)], (3) 
which, calculated for a finite system of linear size L in a pseudo-critical point, scales as 

/3(A„L)~L-V^ (4) 

Usually, in order to see this scaling behavior, one takes into account a sequence of 
values calculated at the pseudo-critical points for some values of Lj, such 
that Li = Lj_i + 1. Expanding 

~ (1 + ^ 1 _ il + ... (5) 

U Li 

one finds a linear behavior of (1 — (3{Li) / (3{Li_i) vs. 1/L for all i. This is the 'linear' 
approximation. Note that the error is 0(1/L^). However, if one can not find a sequence of 
finite systems fulfilling Lj = + 1 (in our case Lj = Lj„i + \/2) this approximation is 
rather crude because of the order of the error. This may be improved in the following way. 
First, let us note that in the following expansion, for small x 

the term x^ is absent and the error is 0{x^). Second, let us put 

X = I^LZhzl^ (7) 
Li + Li_i 

and expand 



Consequently on has a linear dependence of ln[/5(Lj)//5(Lj_i)] vs 1/L: 




(8) 



ln[/5(L,_i)//3(L,)]~^i + ..., 



(9) 



where Lj = (Lj + Lj_i)/2. The main advantage of this approach is a small finite size error, 
what in consequence enables one to examine smaller systems. However, there remains yet 
another problem. It is possible to consider scaling from two pseudo-critical points: (8-18) 
where a; = 1/5 and (18-32) — where x = 1/7. In the third point the expansion |] gives 
rather large finite size correction (x = 1/3) and this point has to be excluded. Thus, one 
can ask whether it is possible to find a scaling relation and to estimate critical exponent from 
two points only? The answer is yes, but the final error will be larger. Since the expansion 
(|) produces an error O(x^), it seems to be especially well suited to this purpose. To test 
this approach we have extrapolated the 1/v exponent from the data for transverse Ising 
model on square lattice ||^ taking into account only two pseudo-critical points: (16-25) and 
(25-36). The estimate oil/v = 1.586(7) obtained by simple linear extrapolation from these 
two points should be compared to the original one 1/v = 1.591(1) extrapolated from 
four points. 

The data collected in Table 1 and 2 enable one to calculate /3-function from Eq. (|^) and 
consequently 1/v from Eq. (P) for each pseudo-critical point of J — J' model; their values 
for Heisenberg and XY Hamiltonians are listed in Table 3 and the graphical extrapolation is 
shown in Fig 3. The values of critical exponents for the same universality classes obtained 
by other authors are also displayed in Table 3. 

One should note that the error in the present approach is larger than in other finite-size 
scaling analysis, but as was mentioned this is a result of the small number of systems with 
required symmetry, see Fig. 1, and the error is still comparable with an error resulting from 
extensive Monte-Carlo simulation []rT],|T2[. 

To conclude, we have presented the results of the investigation of the critical parameters 
for the quantum phase transitions in the spin-^ J — J' Heisenberg and XY antiferromagnets 
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on square lattice. The obtained values of the correlation length critical exponents strongly 
suggest that these transitions belong to the 3D 0(3) and 3D 0(2) universality classes, 
respectively. 
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FIGURES 

FIG. 1. The J — J' model on square lattice. Finite systems of 8 and 18 spins are marked by 
dashed lines. 

FIG. 2. Extrapolation of the finite-system pseudo-critical points as a function of Open 
squares - Heisenberg model, filled squares - XY model. 

FIG. 3. Extrapolation of the finite-system estimates of the l/v exponent a function of 1/L. 
Open squares - Heisenberg model, filled squares - XY model. 
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TABLES 

TABLE I. Pseudo-critical points Ac calculated for two Heisenberg systems of sizes given in the 
first column. The values of the gap A and its derivative A' at the pseudo-critical point Ac are also 
listed. To find them 11 ED data points equally spaced around Ac were fitted to the polynomial of 
4-th order in A in the region of Ac ± 0.005. 5 digits are exact. 



System size 


Ac 


A 


A' 


8 


3.1166 


4.14495 


1.52350 


lo 




z.ToooO 


l.YlOOO 


8 


2.9251 


3.85980 


1.45296 


32 




1.92990 


1.81304 


18 


2.7648 


2.20467 


1.44663 


32 




1.65350 


1.62565 


extra]^olated 


-16(2) 








TABLE IL Same as 


TABLE I, but for XY spin systems. 




System size 


Ac 


A 


A' 


8 


4.9326 


2.51653 


0.786928 


18 




1.67768 


0.848422 


8 


4.7705 


2.39028 


0.770058 


32 




1.19514 


0.942901 


18 


4.6444 


1.44152 


0.788049 


32 




1.08111 


0. 860504 


extrapolated 


4.56(2) 
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TABLE III. Values of critical exponent \/v calculated for two pseudo-critical points for finite 
Heisenberg and XY spin systems and subsequently extrapolated to infinity. For comparison are 
also included values of the same exponents calculated by other authors. 



System size Heisenberg 0(3) XY, 0(2) 



8-18 


1.9841 


1.6220 


18-32 


1.8323 


1.6014 


extrapolated 


1.44(4) 


1.55(4) 




1.46(8) 111 


1.495(5) H 




1.418(6) m 


1.49(1) il 






1.51(1) m 



0.65 




